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ABSTRACT 

Context. Submillimeter spectroscopic observations of comets provide an important tool for understanding their chemical composition 
and enable a taxonomic classification. 

Aims. We aim to determine the production rates of several parent- and product volatiles and the 12 C/ 13 C isotopic carbon ratio in the 
long-period comet C/2004 Q2 (Machholz), which is likely to originate from the Oort Cloud. 

Methods. The line emission from several molecules in the coma was measured with high signal-to-noise ratio in January 2005 at 
heliocentric distance of 1.2 AU by means of high-resolution spectroscopic observations using the Submillimeter Telescope (SMT) at 
the Arizona Radio Observatory (ARO). 

Results. We have obtained production rates of several volatiles (CH 3 OH, HCN, H 13 CN, HNC, H 2 CO, CO, and CS) by com- 
paring the observed and simulated line-integrated intensities. We calculated the synthetic profiles using a radiative transfer code 
that includes collisions between neutrals and electrons, and the effects of radiative pumping of the fundamental vibrational lev- 
els by solar infrared radiation. Furthermore, multiline observations of the CH3OH J = 7-6 series allow us to estimate the rota- 
tional temperature using the rotation diagram technique. We find that the CH3OH population distribution of the levels sampled 
by these lines can be described by a rotational temperature of 40 ± 3 K. Derived mixing ratios relative to hydrogen cyanide are 
C0/CH30H/H2C0/CS/HNC/H 13 CN/HCN = 30.9/24.6/4.8/0.57/0.031/0.013/1 assuming a pointing offset of 8" due to the uncertain 
ephemeris at the time of the observations and the telescope pointing error. 

Conclusions. The measured relative molecular abundances in C/2004 Q2 (Machholz) are between low- to typical values of those 
obtained in Oort Cloud comets, suggesting that it has visited the inner solar system previously and undergone thermal processing. The 
HNC/HCN abundance ratio of ~3.1% is comparable to that found in other comets, accounting for the dependence on the heliocentric 
distance, and could possibly be explained by ion-molecule chemical processes in the low-temperature atmosphere. From a tentative 
H 13 CN detection, the measured value of 97 ± 30 for the H 12 CN/H 13 CN isotopologue pair is consistent with a telluric value. The 
outgassing variability observed in the HCN production rates over a period of two hours is consistent with the rotation of the nucleus 
derived using different observational techniques. 

Key words, comets: individual: C/2004 Q2 (Machholz) - molecular processes - radiative transfer - techniques: spectroscopic - 
submillimeter: planetary systems 


1. Introduction 

Comets spend most of their lifetime in the outer solar system 
and therefore have not undergone much thermal processing. Line 
emission from cometary atmospheres at submillimeter and radio 
wavelengths is a very useful tool for studying their physical and 
chemical conditions and relation with other bodies in the solar 


* Based on observations carried out with the 10-m Submillimeter 
Telescope at the Arizona Radio Observatory, Steward Observatory, 
Mount Graham, Arizona, USA. 
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system (Biver et al. 2002; Bockelee-Morvan et al. 2004). The 
coma structure and expansion velocity can be derived by fitting 
the observed line shapes using a molecular excitation code. In 
addition, mixing ratios of volatiles such as CH 3 OH, CO and CS 
can be compared with observed chemical abundances in proto- 
planetary disks to improve our understanding of planet forma- 
tion processes. 

The composition of comets has been investigated in the 
last two decades to develop a classification based on abun- 
dances of primary chemical species that displays a great com- 
positional diversity (A’Hearn et al. 1995; Mumma et al. 2003; 
Bockelee-Morvan et al. 2004). More than 20 parent volatile 


Article published by EDP Sciences 


A2, page 1 of 13 


A&A 545, A2 (2012) 


species that release directly from ices in the nucleus, in addi- 
tion to several radicals and ions formed by photodissociation in 
the coma, have been detected via ground-based spectroscopic 
surveys at infrared and submillimeter wavelengths and in situ 
measurements. The composition of some cometary ices show 
strong evidence of processing in the solar nebula and can pro- 
vide clues about their place of formation and subsequent evolu- 
tion. The abundance of HCN relative to water has been observed 
to be roughly constant with a value of 0.1% in several comets 
for a wide range of heliocentric distances (Biver et al. 2002). 
Other species show a wide spread of production rates. For in- 
stance, CH 3 OH has been found to have a variable abundance 
relative to water, ranging between less than 0.15% to 6% at 
different heliocentric distances (Bockelee-Morvan et al. 2004). 
There is no evident correlation between the observed relative 
abundances and the dynamical class of the comets. Water is the 
most dominant volatile species and is typically used to determine 
relative abundances. Although water is not directly accessible 
from the ground at submillimeter wavelengths, it has been ob- 
served from space using the Submillimeter Wave Astronomical 
Satellite (SWAS), Odin and Herschel satellites (Neufeld et al. 
2000; Lecacheux et al. 2003; Hartogh et al. 2009) or inferred by 
observations of the hydroxyl (OH) radical at radio wavelengths 
(Crovisier et al. 2002). Herschel ’ s Heterodyne Instrument for 
the Far Infrared (HIFI) is able to determine water production 
rates accurately (Hartogh et al. 2010). Direct measurements of 
water can also be performed using ground-based telescopes by 
observations of non-resonance fluorescence emission at infrared 
wavelengths (Mumma et al. 1986; Bonev et al. 2009). 

Comet C/2004 Q2 (Machholz) was discovered on 27 August 
2004 by Donald E. Machholz (Machholz et al. 2004). The comet 
passed perihelion on 24 January 2005 at a heliocentric distance 
of r h = 1.205 AU and geocentric distance of A = 0.435 AU. 
During its closest approach to Earth at a distance of 0.35 AU on 
6 January 2005 it reached a naked-eye visual magnitude of m v ~ 
3.5 as reported in the International Comet Quarterly. Owing 
to its favorable viewing geometry in the northern hemisphere, 
the comet was extensively observed from the ground at various 
wavelengths. C/2004 Q2 (Machholz) is a long-period highly ec- 
centric comet whose origin is most likely the Oort Cloud ac- 
cording to the comet classification scheme by Levison (1996), 
with an approximate initial orbital period of 118 000 years and 
eccentricity of 0.9994658 before the comet was perturbed grav- 
itationally in the inner solar system (Marsden 2004; Nakano 
2006). These values are close to the perihelion osculating ele- 
ments (see Marsden 2005, JPL Small-Body Database 1 ). Its dy- 
namical classification is still a matter of debate because of the 
strong non-gravitational forces that make a backward orbital in- 
tegration very unreliable. 

In this paper we present high-resolution spectroscopic obser- 
vations of several volatiles from comet C/2004 Q2 (Machholz) 
acquired at the Submillimeter Telescope (SMT). Seven species 
are detected, namely CH 3 OH, HCN, HNC, H 2 CO, CO, CS, 
and a marginal detection of H I3 CN. The comet was observed 
shortly pre-perihelion in January 2005 when it was at a distance 
~0.36 AU from Earth. These observations provide information 
about the outgassing of several molecules and an isotopologue 
of HCN relative to HCN, which is often used as a proxy for 
water in cometary taxonomies. We calculate the CH3OH rota- 
tional temperature of the ground vibrational level from several 
rotational lines and production rates for the observed molecules 
using a radiative transfer code to fit the observed line intensities. 


1 http : //ssd. jpl .nasa . gov/sbdb . cgi?sstr=C/2004+Q2 


Section 2 presents our SMT observations of comet C/2004 Q2 
(Machholz) and the reduction method. In Sect. 3 the radiative 
transfer models and analysis of the observations are described. 
Finally, we discuss the obtained results in Sect. 4. 

2. Observations 

A spectral line survey of primary volatile species in comet 
C/2004 Q2 (Machholz) was made using the SMT telescope lo- 
cated at the Mount Graham International Observatory (MGIO), 
a division of Steward Observatory on Mount Graham, Arizona 
(Baars & Martin 1996; Baars et al. 1999). The SMT has a 
parabolic 10-m primary dish and a hyperbolic secondary reflec- 
tor. The observations were performed with the 0.8 mm double- 
sideband receiver using various spectrometers (Villanueva et al. 
2005; de Val-Borro et al. 2011). Detected emission lines may 
have a different gain response depending upon which side- 
band the lines were observed in. Thus, the sideband gain ratio 
(Gusb/Glsb) deviates from unity and introduces an additional 
uncertainty of ~10% in the absolute brightness temperature cal- 
ibration. A typical system temperature of 150 K was attained 
during the observations. We observed the comet simultaneously 
with the chirp transform spectrometer (CTS) with a bandwidth 
of 215 MHz, and the acousto-optical spectrometers (AOS A, 
AOSB and AOSC) with total bandwidths of 1 GHz, 970 MHz 
and 250 MHz, respectively. The spectral resolution provided 
was 40 kHz for the CTS, and 934, 913 and 250 kHz for the 
AOSA, AOSB and AOSC. One of the purposes of the obser- 
vations was to test the performance of the newly installed high- 
resolution CTS built at the Max Planck Institute for Solar System 
Research (Hartogh & Hartmann 1990; Villanueva & Hartogh 
2004; Villanueva et al. 2006). High spectral resolution is cru- 
cial for resolving the shape of rotational lines in comets and 
study the gas velocity and asymmetries related to non-isotropic 
outgassing. 

Comet C/2004 Q2 (Machholz) was observed near perihelion 
during six nights in the period 13-18 UT January 2005. Here we 
focus on the observations acquired on the first four nights with 
good observing conditions. The data were taken using the stan- 
dard position switching observing mode, where a reference sky 
position separated from the comet by (T5 was observed for the 
same amount of time and subtracted from the on-source observa- 
tions. We used a 30-s exposure for the source and reference po- 
sitions with 8-min scans. A summary of the observing log with 
the total integration times for each line is shown in Table 1. 

The main beam equivalent brightness temperature scale, 
r mB , was corrected for the beam efficiency of the telescope es- 
timated from observations of Mars and Saturn in the interval 
fall 2007-spring 2008 2 and calibrated using the chopper-wheel 
method (Ulich & Haas 1976). The frequency scale was con- 
verted into Doppler velocities in the nucleus frame and corrected 
for the relative motion of the comet with respect to the observer. 
We used orbital elements provided by the JPL HORIZONS sys- 
tem 3 to track the comet and the ephemeris to calculate the po- 
sition and relative motion of the comet with respect to the tele- 
scope. Our observing strategy included pointing and calibration 
observations approximately every hour. The pointing of the tele- 
scope was checked and corrected by dedicated reference ob- 
servations of bright sources like Saturn, the protobinary sys- 
tem W3(OH) and the massive star-forming region DR21(OH) 

2 http://kpl2m.as . arizona.edu/smt_docs/smt_beam_eff. 
htm 

3 http://ssd.jpl.nasa.gov/7horizons 
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Table 1 . Log of the SMT observations of comet C/2004 Q2 (Machholz) in January 2005. 
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(r h ) d 

<A> e 

W 

Beam size 

(UT) 




(min) 

(AU) 

(AU) 



12.96 J 

f HCN 

h 2 co 

4-3 

4-3 

13 

107.2 

1.220 

0.361 

42? 5 

20793 
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22 
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21765 


Notes. In some receiver tunings several lines were observed simultaneously in the lower and upper sidebands. (o) Start times in universal time (UT) 
fractional day values. (6) Number of individual scans. w Total integration time. Heliocentric distance. (e> Geocentric distance. w Solar phase 
angle (Sun-C/2004 Q2 (Machholz)-Earth). 


because they were close to the comet. Typical pointing errors 
were ~4", i.e., about a quarter of the half-power beam width 
(HPBW) at the observed frequencies. 

These observations were very challenging owing to the rela- 
tively uncertain ephemeris of C/2004 Q2 (Machholz) since they 
were carried out about five months after the comet was discov- 
ered, and particularly because of the strong non-gravitational 
forces needed to reproduce the observations afterwards. Its or- 
bit changed significantly during its passage through the inner 
solar system due to planetary perturbations from an initial or- 
bital parameter of a q 1 = 0.000404 AU' 1 to a future value of 
0.001856 AU -1 (Marsden 2004). At the time of the observa- 
tions, the comet was moving with a fast apparent motion of 
~300" per hour in the sky relative to the background stars. This 
introduced an additional term of uncertainty in the pointing ac- 
curacy. Comparing the ephemeris actually used for the tracking 
with the HORIZONS ephemeris calculated from the after-the- 
fact orbit including non-gravitational accelerations, we obtain 
an average pointing offset of 4" . Including the effect of the tele- 
scope pointing error, the offset was at maximum 8". This value 
was considered in our computation of the production rates for 
the observed species shown in Table 2. 

The data analysis was performed using the CLASS analy- 
sis software, which is part of GILDAS package 4 , in conjunction 
with the NumPy and SciPy libraries of high-level mathematical 
tools (Jones et al. 2001; Oliphant 2007). A standing wave ap- 
pears as a baseline ripple in some of the spectra. We determined 
the baseline by fitting a polynomial to the emission-free back- 
ground and subtracted it from the original spectrum. After re- 
moval of the standing wave, the individual scans were averaged 
to increase the signal-to-noise ratio (S/N). Line intensities were 
calculated from the weighted averages of the spectra where the 
statistical weights are the inverse square of the root mean square 
(rms) noise of each individual spectrum. 

The comet’s heliocentric distance was in the range r h = 
1.213-1.220 AU, and the geocentric distance was A = 0.361- 
0.376 AU during the observing period (see Table 1). The HPBW 
of the telescope varied between 207 5 and 22" at different fre- 
quencies, corresponding to 6600 and 6900 km projected on the 
comet, so that molecules from the outer coma contribute to the 


4 http : //www . ir am . fr/IRAMFR/GILDAS 


detected emission. The line excitation in this region is domi- 
nated by collisions between neutrals and electrons, and infrared 
fluorescence by solar radiation. During the observing run the so- 
lar phase angle (Sun-C/2004 Q2 (Machholz)-Earth) ranged from 
42? 5 to 45? 0. We detected seven species (CH 3 OH, HCN, H 13 CN, 
HNC, H 2 CO, CO, and CS), whose spectra are shown in Figs. 1, 
2 and 4-9. The observations were focused on the study of the 
hydrogen cyanide chemistry with the detection of HCN, HNC 
and a tentative detection of H I3 CN to determine the isotopic 
l2 C/ 13 C ratio, as well as other volatile species that sublimate di- 
rectly from the nucleus into the coma and abundant daughter 
species. 

In Table 2 we show a comparison of the mixing ratios of 
the detected species with respect to H 2 O with those derived 
from infrared observations on similar dates (Bonev et al. 2009; 
Kobayashi & Kawakita 2009). For species that are detected in 
several nights during our observing run, the mixing ratios refer 
to the weighted average of the observed spectra. Table 3 shows 
the molecules and transitions detected in our survey with 1-cr un- 
certainties. Line frequencies were obtained from the latest online 
edition of the JPL Molecular Spectroscopy Catalog (Pickett et al. 
1998). Most transitions were observed with the higher resolu- 
tion CTS, which resolves the line shape, and the AOS backends 
simultaneously. Integrated line intensities shown in Table 3 were 
obtained integrating over velocities in the interval [-2, 2] km s -1 
unless another line overlaps within that region. Doppler shifts in 
the comet rest frame were obtained as the first moment of the ve- 
locity over the same interval (£ ; T mBi Vi/ T m B; where v is the 

Doppler velocity and the index i refers to the channel number). 


3. Results 

3.1. Radiative transfer modeling 

We adopted a molecular excitation and radiative transfer model 
based on the Sobolev escape probability method, which includes 
collisional effects and infrared fluorescence by solar radiation to 
derive the production rates (Bockelee-Morvan 1987; Biver 1997; 
Biver et al. 1999). In the outer coma, solar infrared pumping of 
vibrational bands followed by instantaneous spontaneous decay 
establishes a fluorescence equilibrium. To test the effect of the 
infrared excitation, an independent molecular excitation model 


A2, page 3 of 13 


A&A 545, A2 (2012) 


Table 2. Comparison of production rates of detected species relative to HCN and H 2 0 in comet C/2004 Q2 (Machholz) with statistical 
uncertainties. 


Molecule 

2“ 

2/ 2hcn 

2/ 2h 2 o 6 

2/2h 2 o c 

QIQ^o d 


(molec. s -1 ) 


(%) 

(%) 

(%) 


HCN 

(2.26 ± 0.02) x 10 26 

1 

0.084 ± 0.001 

0 

Ul 

1 + 
O O 

0 © 

s> h- 

0.16 ±0.01 

CO 

(7.0 ± 0.6) x 10 27 

30.9 ± 2.6 

2.6 ± 0.2 

5.07 ±0.51 


CH3OH 

(5.5 ± 0.6) x 10 27 

24.6 ± 2.5 

2.1 + 0.2 | 

2.14 ±0.12 
1.54 ± 0.07 c 

1.2 ±0.1 
1.65 ±0.09' 

h 2 co J 

\ (3.87 ± 0.20) x 10 26 

(1.09 ± 0.06) x 10 27 - 7 

1.7 + 0.1 
4.8 ±0.3 

0.14 ±0.01 
0.41 ± 0.02 

0.11 ±0.03 

0.18 ±0.01 

CS J 

\ (1.15 ± 0.04) x 10 26 

(5.1 ± 1.7) x 10“' 

(4.3 ± 0.1) x 10“ 2 



(1.29 ± 0.04) x 10 26 - 7 

(5.7 ± 1.9) x 10- 1 

(4.8 ± 0.2) x 10“ 2 



HNC 

(7 ± 2) x 10 24 

(3.1 ±0.9) x 10“ 2 

(2.6 ± 0.7) x 10- 3 



h 13 cn 

(3+ l)x 10 24 

(1.3 + 0.5) x 10“ 2 

(1.1 ± 0.4) x 10“ 3 




Notes. (<I> Weighted mean of production rates measured at the SMT on 13-16 January 2005 derived from excitation and radiative transfer models 
assuming a pointing offset of 8". (i>) A water production rate of Qh 2 o = 2.7 x 10 29 molec. s -1 was used to derive the mixing ratios, intermediate 
between the measurement on 20 January 2005 by the Odin satellite (Biver et al. 2007b) and the value derived from observations at infrared 
wavelengths on 19 January 2005 (Bonev et al. 2009). (c) Weighted mean of mixing ratios obtained on 28-29 November 2004 and 19 January 2005 
with NIRSPEC on the Keck II 10-m telescope (Bonev et al. 2009). (<7) Mixing ratios obtained on 30 January 2005 with NIRSPEC on the Keck II 
10-m telescope (Kobayashi & Kawakita 2009). 69 Revised CH3OH mixing ratios from Villanueva et al. (2012). 771 Production rates for H 2 CO and 
CS assume the presence of an extended source described in the text. 


based on the publicly available accelerated Monte Carlo radia- 
tive transfer code ratran was used to calculate the populations 
of the rotational levels as a function of the distance from the nu- 
cleus for various molecules and line emission in the cometary 
coma (Hogerheijde & van der Tak 2000) 5 . This code includes 
collisional effects with water molecules and electrons, but ne- 
glects the pumping by solar infrared radiation from the ground- 
state vibrational level. We used the one-dimensional spherically 
symmetric version of the code following the description outlined 
in Bensch & Bergin (2004), which has been extensively tested 
and used to interpret Herschel cometary observations (see e.g. 
Hartogh et al. 2010; de Val-Borro et al. 2010; Hartogh et al. 

201 1)7 

The radial gas density profiles for parent molecules H 2 0, 
CH 3 OH, HCN? H 13 CN, HNC, and CO were obtained using 
the standard Haser spherically symmetric distribution for parent 
volatiles (Haser 1957): 


"pO) 


Q / r Py\ 

4/rr-u exp \ u eX p / 


( 1 ) 


where Q is the total production rate in molecules s -1 . The pho- 
todissociation rate /J p takes into account the dissociation and ion- 
ization of molecules by the solar UV radiation, u e xp is the ex- 
pansion velocity and r is the nucleocentric distance. However, 
second-generation molecules such as HNC, CS and H 2 CO are 
not correctly described by this profile. In the last two cases, the 
density profile of daughter species originating from the photodis- 
sociation of an extended parent source in the coma is given by 
(Combi et al. 2004) 


«d (r) = 


Q Pv 
4flT 2 y eX p/f d — jS p 



( 2 ) 


where /? p and /3 d are the parent and daughter photodissociation 
rates. We obtained the photodissociative lifetimes for CH 3 OH, 
HCN, HNC, H 2 CO and CO from Crovisier (1994) assuming the 
quiet-Sun reference spectrum, and the CS radical lifetime was 

5 The source code is available from ratran ’ s website at http : / /www . 
sron.rug.nl/~vdtak/ratran/frames.html 


taken from Biver et al. (2011). The photodissociation rates were 
scaled by the heliocentric distance of C/2004 Q2 (Machholz) 
at the time of the observations. Deviations from a spherically 
symmetric Haser distribution due to anisotropic outgassing were 
not considered. 

Neutrals in cometary atmospheres are excited by collisions 
with other molecules and electrons, which are the dominant ef- 
fects in the inner coma, and radiative pumping of the fundamen- 
tal vibrational levels by the solar infrared flux. Infrared pumping 
of vibrational bands by solar radiation contributes to the excita- 
tion in the outer coma where the gas and electron densities are 
low (Bockelee-Morvan 1987). Most of the detected emission in 
comet C/2004 Q2 (Machholz) originates from a region in which 
molecules are in an excitation state intermediate between fluo- 
rescence and collision-dominated equilibrium. 

Collision rates with water and electrons were obtained for 
several molecules from the current version of the Leiden Atomic 
and Molecular Database 6 (LAMDA; Schoier et al. 2005). When 
collision rates with water were not available, we scaled the 
collision rates with molecular hydrogen by the ratio of their 
molecular weights. However, the calculated line intensities of 
the rotational transitions in the ground-based vibrational level 
are weakly dependent on the collision cross-sections - of about 
a few percent. 

We used a water production rate of <2h 2 o = 2.7 x 

10 29 molec. s -1 , close to the Odin satellite measurement on 
20 January 2005 of (2.64 ± 0.08) x 10 29 molec. s -1 derived from 
the observation of the fundamental transition at 557 GHz (Biver 
et al. 2007b) at r h = 1.208 AU, and to the water production 
rates of (2.73 ± 0.07) x 10 29 molec. s“' and (2.76 ± 0.08) x 
10 29 molec. s -1 derived from observations at infrared wave- 
lengths on 19 January 2005 (Bonev et al. 2006, 2007, 2009). 
Excitation parameters in the model are the neutral gas kinetic 
temperature, which controls the molecular excitation in the colli- 
sional region, and the electron density. We assumed a gas kinetic 
temperature of 60 K and the rotational temperature derived from 
the relative line intensities of the J = 7-6 transitions for CH 3 OH. 
The electron density and temperature profiles from Biver (1997) 


6 http://www.strw.leidenuniv.nl/~moldata/ 
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Table 3. Rotational emission lines in comet C/2004 Q2 (Machholz) observed by SMT on 13-16 January 2005. 
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Fig. 1 . CH3OH averaged spectrum obtained on 13.95 UT January 2005. The lower and upper x-axis scales represent the frequency of the lower and 
upper sidebands, respectively. Labels indicate detected CH3OH spectral lines listed in Table 3. There is a blend of three A + and A~ emission lines 
at 338.513 GHz (label 5). The observed line at 338.722 GHz is a blend of two E components (label 13). The A~ line at 341.416 GHz (label 14) 
is in the upper sideband close to the A + transition at 338.409 GHz (label 3). A marginal detection of the E line at 338.530 GHz is indicated by 
label 6. The unidentified emission features labeled by (?) are discussed in Sect. 3.2. 


were assumed. Since the electron density in the coma is not 
a fully known quantity, an electron density scaling factor of 
x„ e = 0.5 with respect to the reference profile deduced from 
in situ measurements of comet 1 P/Halley was used. The expan- 
sion velocity is assumed to be constant in the coma, with a value 
of 0.75 kms -1 , obtained from the width of the HCN line. This 
value is consistent with the shape of several lines observed with 
the Institut de Radioastronomie Millimetrique (IRAM) 30-m ra- 
dio telescope at about the same time (Biver et al., in prep.). 

The radiative transfer equation was solved by integrating 
along various lines of sight through the coma that covered 
2.5 x HPBW at a given offset from the nucleus, to obtain the 
brightness distribution in the plane of the sky once the level pop- 
ulations were calculated. The beam-averaged emission was then 
calculated at the distance of the comet. 

Production rates were calculated comparing observed line- 
integrated intensities with the integrated synthetic lines obtained 
with the Sobolev escape probability and ratran codes by fit- 
ting the production rate to the observed line intensity to de- 
rive the relative molecular abundances in the nucleus (Bensch & 
Bergin 2004). Opacity effects are negligible for all lines, there- 
fore the difference in the derived production rates in our mod- 
els arise only from including the pumping of excited vibrational 
bands by infrared solar radiation as an excitation mechanism 
(Bockelee-Morvan & Crovisier 1989). Table 2 summarizes the 
mixing ratios relative to water and hydrogen cyanide in comet 
C/2004 Q2 (Machholz). We show the derived production rates 
with statistical uncertainties for the observed molecules using 
parent- and daughter density distributions in Table 3. Production 
rates are affected by pointing errors of up to 8 " since the 
ephemeris used at the time of the observations was relatively 


uncertain. Therefore we considered this pointing offset to com- 
pute the production rates. 

3.2. CH 3 OH 

Methanol has been observed in several comets with a wide range 
of abundances since its initial discovery in comet C/1989 XI 
(Austin) by Bockelee-Morvan et al. (1990). Methanol rotational 
lines often appear in multiplets at millimeter and sub-millimeter 
wavelengths, allowing the estimation of the temperature and ex- 
citation conditions in the coma. Rotational states are ordered in 
non-degenerate A levels and degenerate E levels corresponding 
to different symmetry states. Radiative transitions between A 
and E species are strictly forbidden, and collision-induced tran- 
sitions are highly improbable. We detected twelve A + , A~ and 
E-methanol emission lines of the J = 7-6 series in C/2004 Q2 
(Machholz) at frequencies listed in Table 3. In addition, there is 
a blend of three emission lines at 338.513 GHz and a blend of 
two lines at 338.722 GHz. Figure 1 shows the observed CH 3 OH 
rotational spectrum between 338.10 and 338.75 GHz with labels 
indicating the transitions listed in Table 3. Several unidentified 
emission features are found at about 338.178, 338.232, 338.277, 
and 338.651 GHz in the lower sideband and 341.169, 341.543, 
341.588, and 341.642 GHz in the upper sideband, which are la- 
beled by (?) in Fig. 1 . The line widths of these features are com- 
parable with the detected CH 3 OH lines and thus are consistent 
with a cometary origin, although no plausible lines were found 
at those frequencies nor the corresponding image band frequen- 
cies in the JPL Molecular Spectroscopy Catalog (Pickett et al. 
1998) and the Cologne Database for Molecular Spectroscopy 
(Muller et al. 2001, 2005). Integrated intensities and velocity 
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Fig.2. Averaged spectrum of the CH 3 OH (13i-13oA~) line at 
342.730 GHz observed on 13.20 and 16.16 UT January with the CTS. 
The vertical axis is the calibrated main beam brightness temperature 
and the horizontal axis is the Doppler velocity in the comet rest frame. 


shifts of the identified CH 3 OH lines with statistical uncertain- 
ties are shown in Table 3. The observed velocity shifts toward 
the blue wing suggest asymmetric outgassing with preferential 
emission in the direction toward the observer, which should be 
studied with a non-spherically symmetric model. Additionally, 
we detected serendipitiously a high-excitation-level methanol 
transition J = 1 3 1 — 1 3qA at 342.730 GHz in the CS/CO spec- 
tra from 13 and 16 January. Figure 2 shows the averaged CTS 
spectrum from both dates. 


3.2.1. Rotational temperature 

The rotational temperature was calculated from the relative in- 
tensities of the individual CH 3 OH lines between levels with 
quantum numbers J = 7-6 and 13i-13 0 A - using the rotation 
diagram technique, assuming that the population distribution 
of the levels sampled by the emission lines is in local ther- 
modynamical equilibrium (LTE), i.e., described by a Maxwell- 
Boltzmann distribution characterized by a single temperature. 
Then the column density of the upper level within the beam, N n , 
can be expressed as 


N u = N 


gu 

Z(T mt ) 



E a \ 

k B Trot / 


(3) 


where g u is the degeneracy of the upper level, Z denotes the par- 
tition function, which is a function of temperature, T rot is the 
rotational temperature, E u is the energy of upper state, k B repre- 
sents the Boltzmann constant, and N is the total column density 
averaged over the beam. This method is commonly used in stud- 
ies of the interstellar medium. 

A linear fit of \n(N u /g u ) versus E u provides the rotational 
temperature as the inverse of the slope and the averaged total col- 
umn density divided by the partition function from its intercept 
at E u = 0. We included in the analysis the blended transitions at 
338.722 and 338.723 GHz starting from the I 2 and 1 -2 E levels, 
and the 338.541 and 338.543 GHz pair from the 7 3 A + and 73 A” 
levels, by assigning half of the observed intensity to a virtual 
level of intermediate energy (see Bockelee-Morvan et al. 1994b). 
The high-excitation-level line at 342.730 GHz was included in 



Fig. 3. Rotation diagram for CH 3 OH lines in comet C/2004 Q2 
(Machholz) including l-<x uncertainties. The column density of the 
upper level divided by its degeneracy in logarithmic scale is plotted 
against the energy of the upper level. Filled (blue) circles are A + and 
A~-CH 3 OH lines, empty (red) circles denote E-CH 3 OH transitions and 
the empty (red) triangle indicates the 3 -cr upper limit on the J = 74 - 64 E 
transition at 338.530 GHz. The solid line shows the best linear fit with 
a derived rotational temperature of 85 ± 7 K. The dashed line shows the 
best linear fit excluding the J = 13i-13oA~ transition at 342.730 GHz 
(energy of the upper state of ~228 K relative to the ground state) with a 
derived rotational temperature of 40 ± 3 K. 


the rotational diagram analysis. For optically thin conditions, 
which normally apply to cometary lines of most volatile species, 
the observed intensity, f T mB du, is proportional to the column 
density of the upper transition level (Bockelee-Morvan et al. 
1994a): 

where A u i is the Einstein coefficient for spontaneous emission 
and d u i the transition frequency. 

The rotation diagram for multiple CH 3 OH lines that fits 
both A and E-CH 3 OH simultaneously was determined using 
a weighted linear least-squares method where the weights are 
equal to the reciprocal of the variance of each line intensity 
measurement. The standard deviation of the fitted parameters 
is obtained from the diagonal elements of the covariance ma- 
trix. Figure 3 shows the best fit for all observed lines and only 
the J = 7-6 series with statistical uncertainties. In contrast to 
water, no reliable A/E ratio was retrieved from radio lines for 
CH 3 OH (Pardanaud et al. 2007). This ratio is almost insensitive 
to temperature, i.e., it significantly departs from the statistical 
ratio only for very low spin temperatures, because of the small 
energy difference of the A and E levels. Retrievals of the A/E ra- 
tio in CH 3 OH for comet C/2001 A2 LINEAR at infrared wave- 
lengths are consistent with T sp i n > 18 K (Villanueva et al. 2012). 

The measured rotational temperature in C/2004 Q2 
(Machholz) is 85 ± 7 K with 1-cr uncertainty from the linear 
best fit of the rotational diagram that includes all the observed 
methanol lines on 13.20, 13.95 and 16.16 UT January (solid line 
in Fig. 3). The high-excitation line J = 1 3] -1 3<)A is too bright, 
indicating that the LTE assumption is not valid for that tran- 
sition. Weighted least-squares fitting methods are known to be 
very sensitive to points that are substantially farther away from 
the linear relation than expected. For a non-thermal population 
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distribution, a different rotational temperature is anticipated for 
the J = 13i-13oA“ line and thus the LTE model is not applicable 
to all the data. In addition, the observations were not simultane- 
ous and variability in the outgassing rate is likely to be present. 
If the J = 13i-13oA~ line is excluded from the rotation diagram 
analysis, a rotational temperature of 40 ± 3 K is derived (dashed 
line in Fig. 3). 

Rotational temperatures of 76 ± 2 K in HCN and 93 ± 2 K in 
fFO were obtained on 19 January 2005 with the Near InfraRed 
echelle SPECtrograph (NIRSPEC) instrument at the Keck II 
telescope (Bonev et al. 2009). Kawakita & Kobayashi (2009) 
found a rotational temperature of 85 ± 5 K in EEC) and 90^° K 
in CH4 on 30 January 2005. Retrievals of the ortho-to-para ra- 
tios in HoO in C/2004 Q2 (Machholz) at infrared wavelengths 
are consistent at 95% confidence level with r sp j n > 34 K (Bonev 
et al. 2007) and r sp j n > 27 K (Kawakita & Kobayashi 2009). 
From the A/E/F relative abundances in CH4 the lower limits 
are 7’ sp i n > 35-38 K (Bonev et al. 2009) and T sp i n > 27 K 
(Kawakita & Kobayashi 2009). The nuclear spin temperature 
is defined as the rotational temperature retrieved from a given 
spin isomers abundance ratio assuming their respective popu- 
lation distributions are in LTE conditions. These observations 
at infrared wavelengths prove the collision-dominated inner re- 
gion of the coma where the population levels are roughly in ther- 
mal equilibrium, and the rotational temperature is expected to be 
similar to the kinetic temperature only in the collisional region 
(Bockelee-Morvan et al. 1994a). However, the derived rotational 
temperature of CH3OH at 338 GHz is likely to be different from 
the kinetic temperature because of relaxation of the population 
levels toward fluorescence equilibrium. We could not determine 
the rotational temperature for other gas species because no other 
simultaneous line observations were obtained. Therefore a ki- 
netic temperature of 60 K was assumed in the excitation model- 
ing for other species. 


3.2.2. Production rate from the rotational diagram 

We derived a CH3OH beam-averaged column density of 
( N)/Z(T rot ) = (2.22 ± 0.45) x 10 11 cm -2 from the intercept of 
the linear fit of the J = 7-6 transitions in the rotational di- 
agram at E u = 0. Assuming that the population of methanol 
levels is close to a Boltzmann distribution, the derived col- 
umn density is a good approximation of the real distribution 
in the coma. Accordingly, we omitted the transition J = 13i- 
13oA - in deriving the column density for the reasons outlined in 
Sect. 3.2.1. We computed the partition function evaluated at the 
rotation temperature by performing a logarithmic interpolation 
between the values from Villanueva et al. (2012). These parti- 
tion functions were computed using two million CH3OH levels 
for the range T = 1-500 K. The resulting column density is 
( N ) = (8.8 ± 3.2) x 10 13 cm -2 including the 1-cr uncertainty 
in the partition function introduced by T mt . A methanol produc- 
tion rate of <2 ch 3 oh = ( 6.8 ± 2.5) x 10 27 rnolec. s -1 was derived 
from the column density assuming a Haser parent molecule dis- 
tribution with a constant expansion velocity of 0.75 kms -1 and 
a beam-pointing offset of 8" from the nucleus, corresponding to 
a mixing ratio relative to water of ~2.5%. On the other hand, 
from the excitation model described in Sect. 3.1, a production 
rate of 2 ch 3 oh = (5.5 ± 0.6) x 10 27 molec. s -1 was estimated 
from the J = 7-6 transitions assuming that 7ki n = T mt and the 
same pointing offset, which corresponds to a relative abundance 
relative to water of ~2.1%. This mixing ratio represents an inter- 
mediate value in comparison with other Oort Cloud comets. 



i’ [km s '] 


Fig. 4. HCN (4-3) multiplet at 354.505 GHz observed with the CTS 
on 12.96 UT January with a 6420 s total on-source integration time. 
The vertical axis is the main beam brightness temperature scaled by the 
beam efficiency of SMT and the horizontal axis is the Doppler velocity 
in the rest frame of the nucleus. 


3.3. HCN 

The HCN transitions are the brightest emission lines in 
cometary atmospheres for ground-based submillimeter obser- 
vations (Schloerb et al. 1987; Bockelee-Morvan et al. 1994b; 
Drahus et al. 2010). Its production rate has been found to be 
in a roughly constant ratio with respect to water, therefore it 
is often used to derive relative molecular abundances. Figure 4 
shows the HCN J = 4-3 rotational line at 354.505 GHz ob- 
served by the CTS. The line is detected with a very high S/N 
of ~100 after about two hours of integration. The intensity and 
velocity shift of the observed HCN emission line intensity are 
shown in Table 3. From the width of the line we obtain an esti- 
mate of 0.75 km s -1 for the coma expansion velocity. Expansion 
velocities in the range v eX p = 0. 5-0.8 kms -1 are typical for 
comets at this heliocentric distance and for this gas production 
rate (Tseng et al. 2007). Because HCN has a short rotational 
lifetime, it is important to model its excitation including col- 
lisional and radiative processes, which control the population 
levels in the outer coma. We derive an HCN production rate of 
(2.26 ± 0.02) x 10 26 molec. s -1 assuming the kinetic tempera- 
ture to be equal to 60 K. The water production rate was assumed 
to be <2h 2 o ~ 2.5 x 10 29 molec. s -1 , which corresponds to an 
HCN/H2O mixing ratio of about 0.08%. This value is moder- 
ately lower than the standard ratio of 0.1% observed at radio 
wavelengths for several comets over a wide range of heliocen- 
tric distances (Biver et al. 2002; Bockelee-Morvan et al. 2004). 

3.4. H 13 CN 

Owing to the strong activity of C/2004 Q2 (Machholz) and 
its close approach to Earth during our observational campaign, 
it was possible to search for weak transitions of rare isotopic 
species. The H 13 CN J = 4-3 transition at 345.340 GHz was 
observed simultaneously with the CS and CO transitions on 
14.20 and 16.16 UT January with a total integration time of 
~640 min, and is marginally detected by averaging the obser- 
vations from both days (see Fig. 5). The inferred production 
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Fig. 5. Averaged spectrum of the H 13 CN (4-3) multiplet at 
345.340 GHz observed on 13.20 and 16.16 UT January with the low- 
resolution AOS. The vertical axis is the calibrated main beam brightness 
temperature and the horizontal axis is the Doppler velocity in the comet 
rest frame. This line is outside the CTS and AOSC frequency ranges. 

rate is (3 ± 1) x 10 24 molec. s _1 . The derived 12 C/ 13 C iso- 
topic ratio of 97 ± 30 is consistent with the standard solar 
value of 89 observed in various comets. Isotopic fractionation 
is very sensitive to chemical and physical conditions where 
the molecules condensed and therefore provide important clues 
about how cometary material formed in the early solar system. 
If comets originate from different regions in the solar nebula, 
they could display variation in their isotopic composition de- 
pending on the local temperature. Additional measurements of 
the 12 C/ 13 C ratio in HCN have been obtained in comets C/1995 
01 (Hale-Bopp) (Jewitt et al. 1997; Ziurys et al. 1999) and 
17P/Holmes (Bockelee-Morvan et al. 2008), which are consis- 
tent with a solar value. The carbon isotopic composition of sev- 
eral comets of different dynamical families also agrees with a 
solar value from measurements of C 2 and CH (see Jehin et al. 
2009; Bockelee-Morvan 2011, for reviews on the isotopic com- 
position of cometary volatiles). 

3.5. HNC 

Hydrogen isocyanide (HNC) is a metastable isomer of HCN and 
has been detected in several comets with a measured HNC/HCN 
ratio ranging between 0.03-0.3 (Lis et al. 2008). Figure 6 shows 
the HNC spectrum in C/2004 Q2 (Machholz) observed by the 
CTS. A production rate of (7 ± 2) x 10 24 molec. s -1 is derived, 
corresponding to a relative abundance of 3.1% relative to HCN, 
a typical abundance for moderately active comets. The origin 
of cometary HNC is still a debated topic. HNC is an abundant 
species that may be produced by chemical reactions from par- 
ent volatile species in the coma. Models have shown that HNC 
could also be formed in the coma via ion-neutral and isomeriza- 
tion chemical reactions (Irvine et al. 1998). Consequently, deter- 
mining the production rate is important to distinguish between 
single parent volatile and multiple precursors for HNC. Since 
we do not have precise estimates on its distribution, we assumed 
the Haser formula for a parent molecule in our computation of 
the production rate. 
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Fig. 6. Averaged spectrum of the HNC (4-3) multiplet at 362.630 GHz 
observed on 14.17 and 15.14 UT January with the CTS. The vertical 
axis is the main beam brightness temperature and the horizontal axis is 
the Doppler velocity in the comet rest frame. The effective resolution 
after smoothing is 200 ms -1 . 


3.6. H z CO 

Formaldehyde is a common molecule in the interstellar medium 
and has been detected in more than 20 comets at radio wave- 
lengths since the first unequivocal detection at millimeter wave- 
lengths of the 3i2-2n line at 226 GHz in comet C/1989 XI 
(Austin) (Colom et al. 1992). Solid-phase hydrogenation reac- 
tions on the surface of the nucleus could account for the pres- 
ence of formaldehyde and have been confirmed recently by lab- 
oratory experiments (Nagaoka et al. 2005; Hidaka et al. 2009). 
Figure 7 shows the H 2 CO 5is— 4i4 transition at 351.769 GHz in 
comet C/2004 Q2 (Machholz) obtained by the AOS. CTS and 
high-resolution AOS data are not available for this transition. It 
is possible that the frequency scale in this range of the AOS can 
be slightly inaccurate. For the acousto-optical backends, the cor- 
respondence between channel number and frequency is not quite 
linear, and the correction may be important and ill-known on the 
border of the band. This line is detected in the lower sideband of 
the AOS at about 300 MHz from the central frequency with a to- 
tal bandwidth of ~1 GHz. The same H 2 CO line was observed at 
the Caltech Submillimeter Observatory (CSO) on similar dates 
(12.1 and 13.1 January) where the velocity shift toward the blue 
wing was found to be about twice smaller (Biver et al., in prep.). 

H 2 CO has ortho- and para spin isomers that differ in the 
alignment of the nuclear spin of the H atoms (parallel in the or- 
tho state and anti-parallel in the para state). We assumed that 
cometary formaldehyde has an ortho-to-para ratio of 3 accord- 
ing to the nuclear spin statistical weights ratio. H 2 CO is a short- 
lived molecule and its spatial distribution has been found to dif- 
fer from that expected for a parent molecule. An extended source 
of H 2 CO has been identified in several comets (Meier et al. 1993; 
Bockelee-Morvan et al. 1994b). Assuming that H 2 CO originates 
either from the nucleus or from an extended source with a par- 
ent scale length of L p = 8000 km, production rates vary between 
(3.87±0.20)xl0 26 molec. s _1 and (1.09±0.06)x 10 27 molec. s“\ 
respectively, as inferred by a model with an 8" average point- 
ing offset (Biver 1997). The first value represents a lower limit 
because it is expected that direct release from the nucleus un- 
derestimates the real H 2 CO production rate. The corresponding 


A2, page 9 of 13 


A&A 545, A2 (2012) 



v [kms 1 J 


Fig. 7. H 2 CO C 5 15 — 4 14 ) transition at 351.769 GHz observed on 12.96 UT 
January with the lower resolution AOS. The vertical axis is the main 
beam brightness temperature and the horizontal axis is the Doppler ve- 
locity in the comet rest frame. This line is outside the CTS and AOSC 
frequency ranges. The frequency scale in the range close to the border 
of the AOS can be slightly inaccurate. 


Gh,co/Gh 2 o mixing ratios are 0.14% and 0.41% for the parent 
source and extended HjCO source production. These values are 
intermediate in comparison with the observed mixing ratios in 
other Oort Cloud comets (Bockelee-Morvan 201 1). 

3.7. CO 

We show the carbon monoxide J = 3-2 line in Fig. 8. The 
line profile is strongly asymmetric with a peak in the red- 
shifted wing. There is some observational evidence that the CO 
lines may be affected by self-absorption effects in CO-rich and 
very active comets, which have been included in our analysis 
(Bockelee-Morvan et al. 2010). For the observed CO transition 
in comet C/2004 Q2 (Machholz) our model predicts an opti- 
cal depth of r ~ 0.2. Therefore, the observed excess emission 
at redshifted frequencies is attributable to kinematic effects by 
the presence of a jet roughly in the antisolar direction since 
the solar phase angle was small during the observations. An 
anisotropic outgassing radiative transfer model with a rotating 
jet-like structure would provide a better fit to the line shape (see 
e.g. Bockelee-Morvan et al. 2009). Based on an isotropic coma 
distribution and using a kinetic gas temperature of 60 K, we de- 
duce a CO production rate of (7.0 ± 0.6) x 10 27 molec. s -1 . CO 
is thought to be a parent molecule with a small additional con- 
tribution from dissociation of CO 2 and we used the Haser parent 
molecule distribution to derive the CO production for simplicity. 

Oort Cloud comets are expected to be enriched in hyper- 
volatile species such as CO and the abundance ratio of FI 2 CO 
and CO has been observed to be correlated in several comets. 
We find a relative abundance of CO with respect to ffCN of 
about 3 1 and a mixing ratio of ~2.6% relative to FFO in C/2004 
Q2 (Machholz). Observations of the hypervolatile CO in comets 
at infrared and radio frequencies show a wide range of abun- 
dances (Bockelee-Morvan et al. 2004). CO abundances observed 
at radio wavelengths range between ~0.4-20% relative to H 2 O 
(Bockelee-Morvan 2011). Generally, these lines are relatively 
weak because of the small CO dipole, although they are the most 
readily observed cometary lines at large heliocentric distances. 



Fig. 8 . CO (3-2) Uansition at 345.796 GHz observed on 14.94 UT 
January with the CTS. The asymmetric line profile suggests a strong 
anisotropic outgassing with preferential production in the direction op- 
posite to the Sun. 


The CO abundance in eight Oort Cloud comets observed by 
infrared ground-based spectroscopy ranged between ~l-20% 
relative to water (Mumma et al. 2003). Thus the measured 
£?co/Gh 2 o mixing ratio lies toward the low end of the observed 
range in Oort Cloud comets. Depletion of CO would therefore 
suggest that C/2004 Q2 (Machholz) has been injected into the 
inner solar system before and is not a dynamically new comet. 

3.8. CS 

The CS radical is believed to be produced in cometary coma 
from carbon disulfide (CS 2 ) and has been observed in several 
comets in the UV and radio wavelengths (Feldman et al. 2004; 
Bockelee-Morvan et al. 2004). Figure 9 shows the observed CS 
J = 7-6 emission line at 342.883 GHz. The line profile is asym- 
metric with a peak in the redshifted side that resembles the CO 
J = 3-2 line shape shown in Fig. 8. This suggests an anisotropic 
outgassing in the anti-sunward direction. It is expected that the 
line is optically thin and absorption effects in the foreground are 
negligible according to our model. To obtain the density profile 
we used a photodissociation rate of fics = 2.5 x 10 -5 molec. s -1 
at r h = 1 AU derived from spectroscopic observations (Boissier 
et al. 2007; Biver et al. 2011). This photodissociation rate corre- 
sponds to a scale-length of 28 000 km, assuming an expansion 
velocity of 0.75 kms -1 obtained from the width of the HCN 
line. A CS production rate of (1.15 ± 0.04) x 10 26 molec. s -1 
is estimated from a spherically symmetric model with direct re- 
lease from the nucleus. Assuming a distributed source in the 
coma with scale length L p = 650 km, which corresponds to 
production from CS 2 (Boissier et al. 2007), a production rate 
of (1.29 ± 0.04) x 10 26 molec. s -1 is inferred. The derived mix- 
ing ratio with respect to water of ~0.04% is lower than the 
typical value measured in other comets at radio wavelengths 
(Bockelee-Morvan 2011; Crovisier et al. 2009a). 

3.9. Short-term outgassing variability 

Several comets have shown periodic variations in the outgassing 
that are associated with the rotation of the nucleus, combined 


A2, page 10 of 13 


M. de Val-Borro et al.: Submillimetric observations of comet C/2004 Q2 (Machholz) 



v [kms 1 J 


Fig. 9. CS (7-6) emission line at 342.883 GHz observed on 14.94 UT 
January with the CTS. A double-peaked profile is observed. The veloc- 
ity scale in the horizontal axis is in the comet rest frame. 



Fig. 10. HCN production rates in comet C/2004 Q2 (Machholz) as a 
function of time on 12-13 January 2005. Each data point represents an 
8-min observation. 


with a long-term seasonal variation with a peak emission close 
to perihelion. In some cases, the rotation periods of cometary 
nuclei have been determined by measuring this periodic out- 
gassing variability using various techniques (see Samarasinha 
et al. 2004, for a review). Rotation periods have been derived 
from the periodic variability of the HCN production rate in 
comets 9P/Tempel 1 (Biver et al. 2007a), 73P-C/Schwassmann- 
Wachmann 3 (Drahus et al. 2010) and 2P/Encke (Jockers et al. 
201 1). To address the question of short-term periodic variability 
in the production rates of C/2004 Q2 (Machholz), we selected 
the HCN and CS observations from one night, 12-13 January, 
which have sufficient S/N. 

The production rates for each individual scan of the HCN 
observations are shown in Fig. 10. The HCN J = 4-3 line 
was detected in each single 8-min scan with an S/N of ~15. 
We observe an increase of the production rate by roughly 
40% during the 2-h observing interval with a mean value of 


O) 



13.05 13.1 13.15 13.2 13.25 13.3 13.35 

Date [UT] 


Fig. 11. Time evolution of the retrieved CS production rates in comet 
C/2004 Q2 (Machholz) on 13 January 2005. Each data point represents 
an 8-min observation with statistical uncertainties. 


(2.09 ± 0.02) x 10 26 molec. s -1 , although the first scans at the 
beginning of the observations were obtained at lower eleva- 
tion starting at 45° and could be affected by stability problems 
in the instrument or changes in the main beam efficiency. The 
HCN evolution suggests a periodic variation in the activity of 
the comet potentially induced by the rotation of the nucleus due 
to the presence of active regions on the surface. The large am- 
plitude of the variation can be explained by an intrinsic activity 
change as one or more active areas are exposed to solar radiation 
during the rotation cycle. This variation is compatible with the 
rotation periods of 9.1 ± 0.2 h derived by Reyniers et al. (2009) 
from broadband filter optical photometry, 9.1 ± 1.9 h reported 
by Sastri et al. (2005) from studying dust fans visible in R-band 
images, 17.60 ± 0.05 h from narrowband CN images reported by 
Farnham et al. (2007), and 17.8 ± 0.5 h by determining the inner 
coma morphology from photometric variations (Manzini et al. 
2012). However, we are not able to perform a detailed periodic- 
ity analysis of the HCN production curve due to the insufficient 
temporal coverage because the total time span of our observa- 
tions is less than one night. Hence it is not possible to determine 
if the brightness increase is a periodic phenomenon or due to 
non-periodic activity changes. 

Figure 1 1 shows the evolution of the CS outgassing rate on 
13 January during a period of 7.5 h. The CS line is relatively 
strong and has one of the highest peak antenna temperatures 
after the HCN line. CS is detected in most of the individual 
8-min scans with an S/N between 3-6. The comet displayed a 
strong variability with a 40% variation around the mean value. 
It is likely that this brightness variability is not a periodic event. 
Small outbursts or activity changes are expected as comets ap- 
proach the Sun, which could explain the observed variability 
(A’Hearn et al. 2005). It is also possible that the strong fluctu- 
ations in the CS production could be partly explained by point- 
ing uncertainties that affect the derived production rates com- 
bined with instrument instabilities. Observations of the CS line 
on 16 January have a lower S/N and only in a fraction of the in- 
dividual scans the line is detected with at least a 3-cr confidence 
level. 
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4. Discussion 

We have observed several molecular species in comet C/2004 Q2 
(Machholz) when it was close to perihelion in the sub-millimeter 
wavelength range with the SMT during six consecutive nights 
in January 2005. These observations led to the detection of sev- 
eral CH 3 OH, HCN, H 13 CN, HNC, H 2 CO, CO, and CS rotational 
lines. Our main goal was to measure the relative production rates 
of several parent and daughter volatiles in this comet and esti- 
mate the CH3OH rotational temperature. We derive an HCN pro- 
duction rate of (2.26 ± 0.03) x 10 26 molec. s -1 at heliocentric dis- 
tance of 1.2 AU, which corresponds to a Qhcn/(2h 2 o ~ 0.08%, 
using a spherically symmetric radiative transfer numerical code 
that includes collisional effects between neutrals and electrons 
and radiative pumping of the fundamental vibrational levels by 
solar radiation (Bockelee-Morvan & Crovisier 1989). In addi- 
tion to the statistical noise, the precision of the production rates 
is affected by the limited pointing accuracy at the time the obser- 
vations were performed. An 8" pointing offset was included in 
the computation of the production rates, which results in ~20% 
uncertainty in the derived values. Mixing ratios relative to hy- 
drogen cyanide and water are listed in Table 2 for the production 
rates derived from the weighted average value of the observa- 
tions on different dates during our campaign. The asymmetric 
shape of the CO and CS line profiles suggests that there is pref- 
erential outgassing from the anti-sunward side of the nucleus, 
while the HCN and H 2 CO line profiles are fairly symmetric, al- 
though the latter is shifted toward the blue wing perhaps due to 
an instrumental effect. Using the rotational diagram technique, 
we retrieve a cold rotational temperature of 40 ± 3 K for the 
methanol energy levels sampled by the J = 7-6 transitions. 

Retrieved molecular abundances relative to water are com- 
parable to those obtained in Oort Cloud comets (Biver et al. 
2002; Crovisier et al. 2009b; Dello Russo et al. 2011). Our 
Ghcn/Gh 2 o mixing ratio of ~0.084% is slightly lower than 
the typical value of 0.1% observed at radio wavelengths, and 
50% lower than those found by Bonev et al. (2009) from av- 
eraged 28 November 2004 and 19 January 2005 observations 
and by (Kobayashi & Kawakita 2009) from observations on 
30 January 2005 using the same instrument. The observed dif- 
ference between our measurement of the HCN mixing ratio and 
those derived at infrared wavelengths is fairly typical. The mix- 
ing ratio of CH3OH relative to H 2 0 agrees with the infrared 
measurement in comet C/2004 Q2 (Machholz) by Bonev et al. 
(2009) performed on 19 January 2005 at r h = 1.208 AU that is 
closer in time with our observations. Considering the pointing 
uncertainty that introduces an error of ~20%, our measurement 
is consistent within 1-cr with the revised abundances retrieved 
from the infrared observations on 19 and 30 January 2005 by 
Villanueva et al. (2012) using a new line-by-line model for the 
V3 fundamental band of CH3OH. The averaged CO production 
rate measured on 13 and 16 January agrees with that obtained 
by Bonev et al. (2009) on 29 November 2004 at r h = 1.493 AU 
((6.3 ± 0.3) x 10 27 molec. s -1 ) within confidence limits. In con- 
trast, the Qco/Qu 2 o mixing ratio of ~2.6% is almost a factor of 
two lower than that derived by Bonev et al. (2009). For a parent 
molecule distribution, the derived g H2 col Gh 2 o ratio is 0.14%, 
and Gh 2 co/Ghcn is 1.7%, which is more reliable since these 
molecules were observed simultaneously. These values are in- 
termediate between those measured by Bonev et al. (2009) and 
(Kobayashi & Kawakita 2009). On the other hand, the values 
inferred from a daughter molecule extended source distribution 
with a scale length of L p = 8000 km are about a factor of 
3 higher. This illustrates that the derived production rates depend 


strongly on the assumed scale length of the parent molecule. 
This is particularly the case for H 2 CO where the beam size is 
slightly smaller than L p , so it may be possible that this value 
overestimates the production rate. An HNC abundance of 3.1% 
with respect to HCN is found assuming direct release from the 
nucleus with a Haser distribution (3.0% if infrared pumping 
of the fundamental vibrational levels is not considered), which 
suggests that this molecule may be destroyed by chemical re- 
actions. This value is compatible with those observed in other 
comets given the dependence on heliocentric distance - except 
73P/Schwassmann-Wachmann and the very active Hale-Bopp 
(C/1995 Ol) (see Lis et al. 2008). 

Outgassing variations induced by the nucleus rotation are ex- 
pected to appear from non-sphericity of the nucleus or the pres- 
ence of active region areas on the surface. The variability ob- 
served in the HCN and CS outgassing rates on 12-13 January 
is affected by inaccurate pointing, which introduces an uncer- 
tainty of about the expected variation caused by the rotation of 
the nucleus. HCN production rates show a uniform brightness in- 
crease over a period of two hours, which could be a periodic phe- 
nomenon consistent with the rotation period of the nucleus de- 
rived using different methods (Sastri et al. 2005; Farnham et al. 
2007; Reyniers et al. 2009; Manzini et al. 2012), but our obser- 
vations do not provide sufficient phase coverage to constrain the 
rotation period. 

Comets are the most pristine objects in the solar system and 
have not undergone substantial thermal processing. Cometary 
ices are more sensitive to thermal processing than dust. Hence 
their abundances provide indications about the formation and 
evolution of material in the early solar nebula. Dynamically new 
comets are expected to be enriched in volatile species, and CO is 
the most volatile component observed in C/2004 Q2 (Machholz). 
However, our observations show a relative depletion of CO and 
and an intermediate-range mixing ratio of H 2 CO and CH3OH 
compared to the range of values measured in other comets (Biver 
et al. 2002). It is believed that the composition of icy material 
on the surface of the nucleus is altered by the exposure to so- 
lar radiation. Thus, these observations provide a hint about its 
thermal/dynamical history and suggest that the comet has visited 
the inner solar system previously. Formation regions of comets 
of Oort Cloud comets vary from 5-30 AU from the protosun 
according to the standard “Nice model” scenario (Morbidelli 
et al. 2008) and can be constrained by observations of the chem- 
ical composition. The relative depletion of several volatiles in 
C/2004 Q2 (Machholz) agrees with the determination of the for- 
mation region in the inner region of the solar nebula compared 
with other Oort Cloud comets, derived from the lower limit of 
the retrieved nuclear spin temperatures and dynamical models 
of the evolution of planetesimals in the solar system (Kobayashi 
& Kawakita 2009; Kawakita & Kobayashi 2009). 

Submillimeter spectroscopy of cometary atmospheres is a 
useful approach for studying the diversity of species that sub- 
limate when a comet approaches the Sun. These observations 
demonstrate the capabilities of the new CTS based on digital 
technology installed at the SMT to support solar system obser- 
vation programs, further established with the cometary observa- 
tions presented in Kuppers et al. (2004); Villanueva & Hartogh 
(2004); Drahus et al. (2010); Paganini et al. (2010); Jockers 
et al. (2011), and observations of the Venusian mesosphere 
in Rengel et al. (2008a, b). Our observations of comet C/2004 
Q2 (Machholz) were analyzed for noise properties, standing 
wave and spectral line properties as compared to model pre- 
dictions. The spectrometer is found to perform according to its 
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original specification in terms of sensitivity, spectral resolution 
and stability. 
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